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ABSTRACT 


The phenomenon of thermally sustained preasure oscil- 
ations wnich frequently occur over liguid helium and 
Miemer Liguiliied gases is modeled and analyzed. The 
model is a distributed system with a linear temperature 
gradient. The interaction between the vavor motion and 
heat transfer is analyzed to derive an expression for the 
time history of the pressure oscillations. 


Cyclic interpretation of the pressure history results 
in a relationship between four parameters which control 
the behavior of the oscillations. The four parameters 
Gemcein information on amplitude of motion, slenderness 
mmeiQ of the tube, characteristic lengths of the tube, 
M@ermel conductivity, svecific heat ratio and viscosity of 
tne gas undergoing the oscillations and boundary layer decay. 


The effect o> changes within the parameters on the 
theoretical behavior of the oscillating system shows good 


Gereement with the behavior of oscillations in previous 
experimental apparatus when similiar changes are made, 


M@aesis Suvervisor: Joseph L. Smith, Jr. 
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Nomenclature 


The following nomenclature will be used throughout 


the text unless otherwise noted. 


symbol 


A 
At 
Ay 


De tewineu ner 
Cross sectional area of tube. 
Surface area for heat transfer. 
Ssuliace area for shear stress, 
Constant coefficient. 
Constant Coefficient. 
Temperature gradient along tube wall. 
Constant. 
specific heat at constant pressure. 
Specific heat at constant volume. 
Viscous drag coefficient. 
Constant coefficient. 
Constant coefficient. 
Diameter of tube, 
Thermal conductivity. 
Length of tube. 


Reduced length for shear stress. 


Length of temperature gradient. 
Length of comoressible volume. 
Mass. 

Non-dimensional perameter. 
Non-dimensional parameter. 


vi 





symbol 


a] | 


Va. 


De teen om 
Non-dimensional parameter. 
Angular frequency. 
Pressure. 

Average pressure, 

Storage vessel pressure. 
Heat transfer rate. 

Gas constant. 

Temperature, 

Cold vapor temperature. 
Warm volume temperature. 
Wall surface temperature. 
Average temperature. 

Times 

Internal energy. 

Velocity. 

Maximum velocity amplitude. 
Volume. 

Average velocity. 

Direction along tube or flat vlate. 


Thullen's maximum vosition of inertial element. 


_ Maximum displacement amvolitude. 


Height above flat plate. 

Maximum heignt above flat plate. 
Thermal diffusivity. 

Dpeciiic Heat ratio. 


Spacial decay rate of momentum boundary layer. 





symool 


mae FY 


seal olay 


We ton 
Kinematic viscosity. 
Viscosity. 
epacial decay rate of thermal toundary layer. 
Devon. 
pnear Strescwon) tubes wali. 
Non-dimensional parameter. 


Non-dimensional parameter, 





I. INTRODUCTICN 
Thermally sustained pressure oscillaticns are 
often found in the vapor trapved within cryogenic equip- 
ment. Particularly suscevtable to this vnhenomenon are 
the tubes and piping which venetrate storaze containers 
of liquified gases. These tubes are characterized by 
meine closed at their outside end which is usually at 
room temperature and open at their other end which is in 
communication with the cold vapor of the liquified gas. 
Bamce it is usually necessary to fill, vent, or measure 
temperature and pressure within storage vessels by means 
Pemee cube Or pine penetration, construction of Cryogenic 
é€quipment without this configuration is nearly impossible. 
Occasionally the oscillations have been used to good 
advantage for such purposes as vapor mixing or liquid 
level measurement. However, in most instances, the 
presence of the pressure oscillations is undesirable. 
When their pressure amplitudes are sufficiently large, 
the oscillations can make vapor pressure He neurones 
erroneous. Since the open end of the tube is in free 
communication with the cold vapor above the liquid, there 
is a heat pumping action to the low temperature region 
maven creatily increaces the boll-off rate of the Jiguified 
gas within the container. 
Encounters with the phenomenon in ervogenic apparatus 
fomvecll documented in cryogenic enzimecring: literature. 
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The phenomenon is most readily observable over liquid 
helium. Investigations by Ditmars and Furukawa (3) and 
Thullen (7) have shown that the pressure oscillations 
may occurmpover liquid hydrogen or liquid nitvrocen, alinouga 
their detection is somewnat more difficult tnan over 
liquid helium. Thullen (7) gives a rather complete survey 
of encounters up to the present time. 

Essentially since the late 1940's, the literature 
bes dealt with the description, detection, and dampins of 
mice .oscililations due to their unwanted nature. Curoslousiy 
enough, it is Lord Rayleigh (5) that is most often quoted 
Mieexplainine the oscillations. In his work ane Theory 
of Sound, Volume II, under the section entitled "Maintenance 
by Heat of Aerial Vibrations," Lord Rayleigh describes 
a form of thermally sustained vibrations encountered by 
glass blowers when blowing a bulb at the end of a long 
narrow tube. Just as found in cryogenic applications, 
mae tube was communicating between temperature extremes 
amaethe tube, itself, certainly must have possessed a 
large temperature gradient along its wall. Qualitatively 
Lord Rayleigh states that, as the gas is moved from the 
cooler regions of the tube to warmer regions, "the adjustment 
of temperature takes time (tralics), and thus the temp- 
Emeture of the air deviates from that of the neighboring 
parts of the tube, inclining towards the temperature of 
that vart of the tube from (italics) which the air has just 


come. From this it follows that at the phase of greatest 





condensation heat is received by the air and at the 

phase of greatest rarefaction heat is siven up from it, 

and thus there is a vendency to maintain the vibrations." 
Lord Rayleicn further points out that the budesandeine 
great range of temverature is necessary for the maintenance 
oretnesvibratiion. Wiithingermyomenticweoparatus Cons temetnon,. 
itew Conductivity meterial such as (¢lasc cure -nicken, 

or stainless steels are used for the tubes. Thus, a 

severe temperature gradient must exist near the cold end 

of the tube; and the warm end of the tube and associated 
valving act as the bulb of Lord Rayleigh's high temperature 
phenomenon. 

Extensive exverimentation has been conducted to 
femermmine the types of materials and geometric confisgurat- 
fons which will foster and sustain the vressure oscillat- 
lons in a room temveratureto cryogenic temperature environ- 
ment. Clement and Gaffney (2), Ditmars and Furukawa (3), 
Bannister (1), and Thullen (7) have in their own ways 
measured frequencies, pressure amplitudes, and devised 
methods of damving the oscillations, From their experiment- 
Jons ani reports it is evident that the tube diameter, 
length, anc magnitude of the longitudinal temperature 
gradient in the tube wall are the most important ohysical 
characteristics of the apparatus. 

Analytical investigations have yet to lead to a 
Satisfactory theoretical explaination of the phenomenon 
associated with cryozenic apparatus. Norton and Mulenhaupt 


(4) have produced e computer program to simulete the 





pressure oscillations. In their vorocram the Continues 
temperature gradient was approximeted by a large number of 
constant temperature regions. The programseeve encoumasune 
results when comvared to exvnerimental measurements of the 
pressure oscillations. Thullen (7) proposed a lumped 
parameter model in which he wes able to form two dimension- 
less variables, one containing frequency, friction and 
temperature gradient information, and the other containing 
amplitude and heat transfer information.. From this 
approach the exact behavior of any distributed system 
Gem not be predicted; but, it does point out the controlling 
variables. Both the approaches made by Norton and Mulen- 
haupt and Thullen have added significantly to the theo- 
retical understanding of thermally sustained vressure 
oscillations in cryogenic equipment. However, due to the 
approaches taken, neither theory was able to capture the 
inherent distributiveness of the system and the continuous 
Hemeuee Of the temperature gradient along the tube wall. 
Pot cmuAcmimteOLiONmel Udicmal Lemme: Cet. e samodcnl 
for the pressure oscillations which encompasses the in- 
portant points of previous experimentation, and analyze 
the model in such a manner as to retain as much distribut- 
iveness of the eeeten and continuity of the temverature 
gradient as possible, while maintaining a manazeable 


theoretical discription of the phenomenon. 





Il. ANALYSIS 


1, The Model. 

Tne model chosen for the theoretical analysis of 
the thermally sustained pressure oscillations closely 
approximates the geometric configuration often described 
in.the literature. Common points in these configurations 
are the closed volume at the outside end of the pene- 
trating tube, the longitudinal temperature gradient along 
the tube wall, and the free communication of the cold 
vapor with the oven end of the tube. The system for 


analysis is illustrated in Figure 1. 


mp, vontrol Volume 
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Liquified Gas 


Figure 1. Model For Analysis 
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It is important that the tenperacurc sa. cdaewe 
alongs the tube wall be of sufficient strength to allow 
heating of the vavor uvon compression and cooling upon 
expansion. The maximum amplitude of the oscillations 
16 limited due to viscous, friction. Thus the mode imo: 
meal systems has all the ellements Vo wlosuereencus coy msome 
and sustain, vressure oscillations. Also, the distributed 
Mewure of the system and the Continuity of tne vemverziuve 
gradient have been preserved, 

Within cryogenic apparatus the oscillations are 
Brontaneous in oriein; and, after a period of =time, their 
amolitudes grow to a maximum value. The analysis of this 
model will assume that a sinusoidally varing steady state 
condition has been reached within the system. The re- 
lationships between time, motion, and heat transfer will 
meenwcombined to define the vressure as a function of time, 
physical characteristics of the vapor, and geometric 


configuration. 


paver ivation of tne Time and Pressure Relationshin. 
When the model is considered as an open systen, 
the oscillations of the vapor within the tube cause 
energy to te added or lost to the system by virtue of 
Hees entering or exiting the open end of the tube, and 
by heat transfered to the vapor because of the temperature 


gradient along the tube wall. The first law of thermo- 


dynamics for an oven system applies to the control volume 





within the closed tube: 
‘e| « 
ae = Cpliae + @ (1) 


The time rate of change of mass entering the system 


across the open énd of the tube is given by; 


dm _ 
an = pAv(t) (2) 
where v(t) is the average velocity of the fluid particles 
at the end of the tube at a given time t. 

The verfect gas relations and the definition of 
anternal energy give the rate of change of pressure in 


terms of the rate of change of internal energy: 


P = URT 
_ Ce) 
i = mCyT : (4) 
_ 
r= Hu (5) 
P= CU | (6) 
where C = RL 
CyV 


Since both pressure and internal energy are functions 


of time: 


_ cau 
oF (7) 


ae 
ctitd 





When equation (7) is substituted into equation (1), 
the rate of change of pressure is given as a function of 
Lime : 


a = G [cpApTyv(t) + a(t)| (8) 


Density is not a variable of primary interest; and, 
over a cycle of the oscillations, its variation is small. 
Therefore, the following simplifing assumption and sub- 
stitution is made in equation (8): 

5 E (9) 
f Bete 
where P is the average pressure at an average temperature, 


T, during a cycle: 


ap _ a [ CpPA,T1 : 
st > o| SEE ent) + 13) (10) 


When v(t) and Q(t) are known, equation (10) can be 
integrated to give the pressure as a function of time. 

Since the usual tube diameters are large in comparison 
to the thermal and momentum voundary layer thicknesses, 
the derivation of the velocity term, v(t), and the heat 
transfer tern, Q(t), in equation (10) will be carried out 
in a rectangular co-ordinate system. ieoere ee A shows the 
rangce of tube diameters for which this avvroach will be 


void. 





3.  Soltittienmot (the Wedioc i aicrm argue 
Schlichting (6) presents the boundary layer solution 
for the flow of a fluid near an oscillating filet plate, 
For the purvose of this analysis, Schlichting'’s solution 
will be modified to fit the boundary condition that the 
fluid is sinusoidally oscillating over a stationary flat 
plate. Accordingly, the velocity distribution in the 
fluid above the flat plate as a function of height above 


the plate, y, and time, t, is given by: 
eet uge” eos (nt - Sy) - Usgcesems (11) 


where § “V8 is the spacial decay rate of the momentum 
boundary layer. 

Bie average velocauy of the Pluidgenterinzstmcecupc 
is obtaired by assuming that the average velocity of the 
fluid over the flat plate to a given distance is analogous 
Gemciat Entering tne tubes 


m 
v(t) = fouty,eyas (12) 
Im 


where yp corresponds to the tube radius, d/2, and is 
sufficiently large so that in equation (12) e?7<e1, 

When the average velocity is determined in this manner, 
ma@estinal Soiution of the velocity term is: 


hy 5 |G ee nt Cos nite 1 sin nt | (13) 
25 Yn ody, 
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4, Solution of (he Heat eeucser Téecnce 

The model assumes that the only heat transfer to and 
from the vapor is in the area of the tube where the 
temperature gracient Cxilecs cs ic mecct Ouse mc Olnlncumons 
the tube is assumed to be adiabatic with wespect acon une 
Memor. AS wish the velocity Venn derivate copmecicet lon 
plate avvroach is assumed to be anplicable. 

To, beeen, anvenergyevalance is made onvauconure: 
volume of unit devth vositioned above the flat plate. 
There is a vertical heat flux through the control volume 


and a horizontal enthalphy flux. Figure 2 illustrates 


mec energy balance, 





-kax82 , 3 ~kdxOl)a Unit Depth 
oy zy oy y Control Volume 
[oa RC ae a 
r udyc_T + 
pudye,? ——> ! oe 
7 | | 3 
i 1 Sx (pudyc pT )ax 
uf | 
~kdxeT 
x oy 
Tw 
a 
E 1 


Figure 2, Energy Balance For Heat Transfer Derivation 
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The energy balance on the control volume yieles the 
folliowins differential equatron: 
oT > 4) 
Pfeyel = eo Fo alae te > (1 
ot “Sy 
Equation (14) may be rearranged into a diffusion equation 


Waeciuaiveloci ty type fone ine unc ikea: 


Labs wot — uly,t)ar (15) 
ZW ot aye ax 


where 3*= mena ol= K_, 
c= vince 

When equation (15) has been solved for the temperature 
Gistribution in the fluid above the flat plate, the 
temperature gradient at the wall can be found by differ- 
entiating the solution with respect to y and setting y 
equal to zero. This will enable the heat transfer term 
to be realized as a function of time. 

The following assumptions are made in reguard to the 
solution of equation (16): 


a. A linear temperature gradient exists in the wall 


so that the wall surface temperature is given by: 


T, = ax +E (16) 
b. The same linear temverature gradient as in the 

wall exists within the gas so that the form 

OL the, lLenperature distro onm mn Lhe was above 


Seo. mimOmes uc Dlabe 1 See) Vemmby . 


Teves) a=) 1 (Vere. eres (17) 





c. There is neglible heat conduction In the eas 

in the x-direction. 

d. The velocity distribution is given by equation 

(11) and the average velocity is given by 
equation (13). 

Equation (17) must also fit the boundary condition 
that et the wall the wall surface temperature and the 
sas temperature must be equal. 

When equation (17) is substituted into equation (15), 
the differential equation which determines the temperature 
distribution above the plate is: 
df yet 


Sp tO ay? 7 BUlysA) (18) 


SS 


where the function f(y,t) 1s subject to the boundary 
Sendition that at y = 0, £(0,t) = EB. 


The solution equation (18) is of the form: 


f(y,t) = ae A,sin (nt - dy) + A,sin (nt - dy) 
+ B sin nt + lee nt 
em D,sin(nt - fv) + D,cos (nt -Zy) 


+ constant (19) 


Appendix B shows the determination of constants associated 
with equatton (19). 
The complete solution to equation (17) giving the 


temperature distribution in the gas above the flat plate is: 





Ty, t)) = ee xe e- Sy sin(nt - dy) 


atdS* on 
ri 
+ u,egsin nt 
n 

“Fy 

~ 2u,add 9/n e?Jsin(nt - Fy) (20) 

| 2tde- n 
Ce 
where f# =1f_n_, spacial decay rate of thermal boundary layer 


Oot 3 
S =\f_n_ , spacial decay rate of momentum boundary 
2v Maver. 
The heat transfered from the plate to the oscillating 
vapor is given by: 
Q=-koT) At (21) 
oyjy =0 
When equation (20) is differentiated with respect to 
y and evaluated at y = O and the result substituted into 


equation (21), the heat transfer term, oF is given by: 


Q = -kad'ugd lots oa! - | jes" nt + cos at | (22) 
2epd*— éf n 
a 


SB. Final Solution of the Pressurs-Time Relationship. 

With the results of sections 3. and 4. combined | 
meen the differential equation of section 2., the final 
derivation of the ges pressure as a function of time can 
be determined. 


Eouation (10) of section 2. is: 





dP _ n~lcuPA 
Eee ee 


2 


waa ih cs Ge) 








Equation (13) of section 3. is: 


v(t) -Seule 








LSS 2 ai coset 4 salon 
20Vn In 


Equation (22) of section 4. is: 


Q = “kad du 2axSfF _ 4||sin nt + cos nt] 
2as* - 2, nD 


When equetions (13) and (22) are substituted into (10) 
and equation (10) is integrated with respect to time, the 


gas pressure as a function of time is: 




















Gt). C CPA 2 ud aa } Sin nt - 1. cos nt 
Ro \y\F [Led ym 25m 
= [sad uf 2A SF Ot } rece nt + sin nt 
2A ff- niin n 
a 
vy kS (23) 


where Po. is the constant of integration or pressure at 
time t=0. 

With the following cefinitions the compressible 
volume vrecsure given by equation (24) can be made non- 


dimensional and more manageable for practical applications: 


Let Nw eel 
€ 0 33 9m 
n 


No = oe d* 
nh 





ye) 


Dae, 
At 
i CoPA 
15 
Rg = kaAd'S a” 
Nn 


Upen the substitution of the above detimiitaenc- 


equation (23) becomes: 


P(t) _ [HE ly IN = I 
| 3 (No - 1) fs = | Sin nt 


CBtu,/ 7 
-jutt, — [Ny - 1 
Bs No 4] CoS ne (24) 


In a likewise manner the average velocity, v(t), can 


be made non-dimensional: 


vit) = (Nj - 1)cos nt + Nosin nt (25) 
0 





Gey ORES Serle 


With the pressure history of the warm volume and 
the motion of the mass taking part in the oscillations 
now Known, a dynamical approach to the system will be 
used to derive the relationship between the varameters 
Woeeh Controi the ostiilations ime smede 1 60 muiem owe com 
for this investigation will be similar to that used by 


Thullen (7) and is shown in Figure 3. 


—2| a fp 
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| 
Cw 
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al), Temoerature Shear Stress 
Erour le on Tube Wall 


Figure 3. Model For Cvelic Interpretation 
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Newton's law of motion 2epplies to the inertial 


element of the models: 
mk + Dk - A[P(t) - PQ(t)] = 0 (26) 


The transient solution to equation (26) is not of 
interest and can be avoided by the cyclic integration 
of equation (26) with resvect to displacement. The 
net work done in accelerating the inertial element through 


eecyclé is zero; 
gf nikax = © | (27) 


The two remaining terms of equation (26) are now 
equated and point out that the cyclic work done on the 
volume to the left of the inertial mass is equal to the 


cyclic work done against friction: 


$ a[P(t) - Pe(t)] ax = g Diax (28) 


The identity ax = a dt = v(t)dt is substituted 
t 
into equaticn (28) so that the cyclic integration can 


be verformed with respect to time instead of displacement. 


 A[P(t) ~ Pelt] vit)at = $Dv2(t)at (29) 


If the exit pressure is constant, Pe(t)-= Bo, the 
following enequality must be satisfied for the oscillations 


to be self sustaining: 


ga[P(e) - Be] v(t )at = @pv?(+)at (30) 





18 


In substituting equations (24) and (25) into (30) 
and verformimsethe cyclic integration tie tel lors 


relationship is establisned: 


ABC (Ny, = 1] >(n. -~ 4)2 N 1 
“nD|No = 4 Za IN = ) = e (31) 








Thullen's (7) non-dimensional parameter containing 
frequency, friction, and temperature gradient information 
is as follows: 

tp = Pope @ ens, (32) 

Weone substitution of the following detinivione 


of the variables in equatiion (32), the Thullen parameter 


pecomes3 
Po = pRTy 
A = Trae 
8s = 1dal' 
b =AT 
Ty 
W = n 


(33) 


oa 
+] 
| 
- ce 
Ld 
Ou 
___ 
ae 
> 
ry 
a 
TT 
ww 
ry js 
| 
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Upon subsitituion of the actusl veriables for B and 
C in the non-dimensional grouping ABC in equation (31) and 
nD 
making the following definitions, this non-dimensional 


group is seen to contain Thullen's parameter plus additional 





information: 
N®¥ = ABC = AkaA'So R_ (34) 
nD n-D CyV 
A = fa? 
UAT 
A' =fTfdl 
=q7d= 1" 
ue 
of =Vin=Vnr=pVa_ 
au 
N¥ =) d2]) pRato}| 1 || _ 4k (35), 
air || 27°" of let 
NY = Npf_4ark (36) 
cyned Cynzd “(up 


iamorcer tO Carryout the scyclic wniccretom on 
equation (29) with a variable exit pressure, Pe(t), 
further definition .of tne exit pressure term and the 
damping coefficient term, D, are required. 
First, the exit pressure is examined: 
b(t) ND Pe (t) 
——=~ + v(t) 
Outflow: P(t) = Po 
Inflow: Pe(t) = Py - Apve(t) 
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Gaia 


$ ae (t)v(t)at = faze(t) 
Cyc Le 


Vv 
nflow $ 


+ fare (t)v(t)at (37) 


out flow 4 cycle 


Tue outflow half-cycle contributes nothing to the 
cyclic integration. In the inflow half-cycle the sign 
of the velocity is negative, Therefore stnuce we vce 
integration of equation (29) reduces to: 


$ a[Pee) ~ Pa,| v(t) at - | aap? (t)at =f v2(t)at (28) 


inflow $ cycl 


Second, the damping coefficient term can be determined 


through a pipe flow analogy involving shear stress$ 
q 
7 Ee 
r= (P, - Po) A =twhy 
tw = du 
FY syly = 0 


The rate of change of velocity with respect to y 


4s found from the velocity equation (11) 


OU ana. 
+h Loe J uo| sin nt +cos nt] (39) 


The damping coefficient is now defined as: 


D= -F = A, Sau [san nt + cos nt| (40) 
y(h) eet | 


Vv pat 


The right half of equation (38) now requires the 


eyelic intezration of the following term: 


ASpu ee nt + cos nt} 2 dt 
; Wy SSR ea v<(t)dt (41) 





ea 


With the above substitutions for the exit pressure 
and damping coefficient terms and subscaucmmeev clic 
integration of equation (38), the relationship between 
the parameters which contribute to the sustawieis Cree 


pressure oscillations can te established: 


3/2 
ABC uSTr 1) 2 Pau Z 2 
20s ioe 2 phué | (No - 1) + NB | = 
1/2 
TTA [(2 Be ee NS (42) 


Upon substitution of the following definitions 


and previous values of B and C, equation (42) tekes 


the form3 
dp -(F he aS oA 
No = 1 of 1 | 
v 


(No _ 1)¢ + ne 


“Ss 
T 


[st - 1° + (zt)? 


Uo = NX_ 


> 
| 


= 7aL', where L' is the reduced length (L - 41') 
of the tube since the shear stress is 
not uniform over the entire length of 


the tube of Figure 3. 


Vi 
Bye 


a, 
T 
he 
Zits 
\ 


7 -1/6 
aE x, = Beas tgs yf 2/7 -¥f (43) 


Wi 





22 


Figure 4, eghows the non-dimensional displecement 


amplitude, : a5 Xo: plotted against the driving force 


parameter, kaTa(o - ae » With lines of constant Us ° 
L'Ligné 


The graphical representation of equation (43) gives 
indication of a given system's performance and how that 
performance will change if the various parameters are 
altered. The manner of change suggested follows closely 
Mauer woich has so often been described in literaturewen 
thermally sustained pressure oscillations. Due to the 
lower limit on tube diameter imposed by the theoretical 
approach taken, the lines of constant U/ in Figure 4, 
mee in agcrement with the theory for Y = 1. Values of 
Us = 1 are shown in dashed lines for completeness. 

Of particular interest is the linear relationship 
between the displacement amplitude and the slenderness 
ratio, L'/d. Bannister (1) noted this in his experimentation 
with oscillations over liquid helium. He concluded that 
the pressure amplitudes are directly proportional to the 
Slenderness ratio of the gas column. The Bannister data 
was used to calculave values of Vf and the eriving Toree 
parameter. The calculated operating points are plotted 
on Figure 4, Apvendix C tabulates the constants and data 
used. Figure 5. shows the calculated amplitude parameter 
plotted against the measured maximum pressure amplitude. 
The actual slenderness ratio of Bannister's tubes is also 


plotted against the maximum pressure amplitude to show 
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the linearity comparison. 


Blan 600 
5Co 
Calaulated Actual 
Amplitude Slenderness 
Parameter 2 400 Ratio 
L/d 


1 a-/ n 
3D" 25 %o 300 


200 


-—--—-- Actual L/d 


Calculated oe 


Amplitude 





0 20 40 60 


Measured Pressure 
Amplitude (ATM x 107) 


Figure 5. Slenderness Ratio Comparison 


Thullen (7) also found a non-dimensional displacement 
amplitude perameter in his enelysis: x, /21. The char- 
acteristic lence was the distance between two ideal 
Meee exchangers in nie model; and, it was viewed as "the 
distance traveled by the ges after a flow reversal before 
Significant heat transfer begins. However, the length 
was not well defined for a distributed system. The 
results of the present analysis indicate that Thullen's 


"1" might have the following value for a distributed system: 


t 
! 


= 
il 
We) 

Quft- 

sl¥ 





25) 
Iv. CONCIUSICNS 


A greater understanding of thermally sustained 
pressure oscillations in liquified gas apparavuc is 
possible througn a study of vapor motion and heat transfer 
interaction. From such a study a time devendent exvuression 
hor tne oressure history of 2 vapor isa ubeec ner 
derived. With cyclic intervretation the relationshin 
anong the controling parameters of the system can be 
established. 

When changes are made in the parameters, the behavior 
Samune theoretical relationships appear to fellow closely 
that which occurs when similiar changes are made in 
experimental apvvaratus. 

The rectangular co-ordinate system from which the 
motion and heat transfer interaction was derived places 
a lower bound on the diameter of the tube for which the 
theory applies. Therefore, an analysis with the same 
Beperal approach to the motion and héat trensier, but in 
Cylinderical co-ordinates, would be of value. It is not 
anticivated that this would alter the basic reletionshinv 
between parameters controling the oscillations; however, 
ae cylinderical system might provide more exact numerical 
constants. Appendix D suggests an approach to the analysis 
of the phenomenon of thermally sustained pressure oscillations 


iia Cylinderical co-ordinate systen. 


























BPPE MN Die A 


In order to show that the tlat plate ssnale= as 
yaletad for the*results desired, it is sutficicutintoucme. 
that the svacial decay rate is small in comparison to 


the tube racius, 


_ on 
J = By 
Let Jy = y/nr 
~ =Ver 
n 
5 x 1074 ft°/sec 
50 rad/sec 


meen A. = 0.0535 in 
e7¥/a 
1.00 
Sie 
050 


025 


12 2a 3BrA4 y 


Pemee, tubes with radii>O.125 im. will be in tne 


menee Lor which the theoretical approach is valid. 


ey 





Toe determination of the constants asscciated with 
equation (19) whica is the solution to the differential 


equation (18) is as follows: 


Equation (18): t af SF 
ake yoy tot, - au(y,t) 


where u(y,t) = ue" *Veos (nt = oy) = Mie ciocm rt 


and f(y,t) must satisfy the boundary condition: f(0,t) = E. 
Equation (19): 


cos(nt - Sy)] 


LSE J) eV [aj sin(nt aot) ey A, 


+ B,sin nt + Bycos ne 
+ e~?¥[D,sin(nt - fy) + Dcos (nt - fy) 
+ constant. 


foe the form of the function f is: fly,t) = f, + fo + Constant 


Seek the solution of equation (18) by fy! 


-§ 
of =e Yl yncos(nt - Sy) - Aonsin(nt - Sy) 

+ B,ncos nt + B,nsin haar 
rf, = 2 f7o~ FY [-4 sin(nt - dy) + A,co3(nt - by)] 
oy 2 | 


Upon substitution of these partial derivatives into 
equation (18) and with like coefficients of the sine end 


cosine terms equeted: 


= U92a 
A, = —~ 
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as. 


2 
An(2eed~- B) = 0, Ne ae 
By = Uoas 
Be =20 


-o 
fl vente wore  sin(nt -dy) + uoatcos nt 
2a Je z, a 


However, £1,aio0ne, will not fit the boundary coud, 
Hence, an additional function must be added to f, and 
that function must be a solution to the homogeneous part 
of equation (18): 1 df _ got 
zot” aye 
fo(y,t) = ef [pysin(nt -¥y) + D,cos(nt - ¥y)| 


+ constant 


eto = ne~*¥[—Dysin(nt = ay) 2 Dycos(nt - fy)| 
2 
oto = 2 Fen ¥|-Dosin(nt -¥y) 12 D,cos(nt - fy) 


Upon the substitution of the partial derivatives 
into the horogenous differential equation and with like 
coefficients ecuated, fp becomes: 

Dol2etf*- 2) = Di(2e F= B) = 0 
If Dj 4 0 and Dy # 0, then (2e¢f*- B) = 0 


Therefore, " Va 


2ot a" 


Q 
ion 


With the boundary condition now applied to f(y,t), 


the coeffients D,;, D, and the "constant" of fo are determined: 





constant 


LHeGen Gio 
constant 


Do 


D, 
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30 
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APPENDIX C 


The following constants for helium gas (-200°F) 


were used in the Bannister data fit of ecauation (43): 


of Daye os 1074 ft©/sec 


yu 
) = 4.0 x 1074 ft2/sec 


Sas 107° lbm/ft-sec 


k = 2.2 x 107° BTU/sec-ft-oFr 
y, 


1.67 


Data from Bannister's exverimentation: 


n (cps) 4 (cm) L (em) P (atm x 10°) AT (°F) 


a 5e0 287 149 56 480 
Ce 5554 ~ 482 149 28 A480 
Oo 33.4 ~605 149 12 480 
A iene 1.240 213 12 480 


Estimated from apparatus geometric configuration: 


L' (cm) 1' (em) 


x 10> 89 
O 105 89 
m 105 89 
oe 105 89 


on 





APPENDIX D 


The cylinderical co-ordinate analysis will proceed 
with analogous steps to those performed in rectangular 
co-ordinates. Figure 5. illustrates the system for analysis, 
Schlichtine’s solution for the velocity dictributionses 


an oscillating flow through a pipe is: 


Navier-Stokes Equation: 


gu _.1_ oP Ou a 
at Ff ox iow + £98) 


Assume that the pressure gradient is caused by a harmonically 


moving piston: 


be 
i hae 


= Koos bt 


Sp fn 
ro 


Use complex notation and attribute physical significance 
only to the real parts; 


- 332 - Kent 


mocume the velocity function has the form: 


u(r,t) = f(rjetnt 





Figure 6. Cylinderical Co-ordinate Model 
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The solution to the Navier-Stokes equation which 


gives the velocity distribution 4001 tiem cee 


u(r,t) = -ikeint I 7 en 
ig 


Jo (RS) 
where 4 =V=22 


The temperature distribution in the vipe is assumed 
to have the following form: 
iG) Saree) amare 


Boundary condition for the temperarure distribution is: 


BRE t) = £ 


Macwaitferential equation resultine from enersy consideravien. 


is: 
1 O2) Rey |e 
23 =a selr3] au(r,t) 
The ultimate differential eauation to be solved is: 
of tar sr 
<= = of = + 3 -~ au(r,t) 


a8 ror 


Assume the solution of the differential eaouation if of the 


Horm : 
ie (haste) memoar Gc te)) foestcylcente) eer. (6 
where 
, int i ine 
f Pde eel eas RE el 
ere? ; | Jo(Rd) 
ae int S J tet) | a . - _ = 4 
Cay Gen iey, = ales | To(a¥) + C, where - Vein 





The solution is: 


| 


| _ _ ak 1 _ Io(rd)] cint 
Bal Far | ace - 


Ser, t) a8 AS A +2 _ |" ~ Jo(r8)| eint Jf 


*(- “St Jo(RF¥) 
The temperature distribution within the pive is: 


T(r,t) ={- | - agizt) _ a . dole 4) Beton: 
Jo(Rd) A Jo(RZ) 


+ E+ ax , where B= 3k , b= 1 SAH, 0 ad? 
-:) 


The heat transfer from the wall surface to the vavor 


in the pipe is; 
ine oT 
Q k Trl 37) = R 


Therefore, the heat transfer term sought is: 


Noe @ f_§o1(Ré) ot int 
Q =-2rk1R A | 35(R3) +07 7027) e 


The average velocity of the vapor entering the pipe 
R 


v(t) = I efru(r,t)dr 
If RE 


mee 2 Jo(Rd int 
eS [ - Be ee] — 


Now both the velocity term and heat transfer terms 


is3 


are functions of time and can be substituted into equation 


Oe) of section II.2,: 


if = o [ends v(t) + To) 





5) 


Upon substitution and integration with respect to time, 
and attaching significance to the real part of the result, 


the pressure as a function of time is: 


P(t) - Ree g ({e08s EIFS) [1 - 93 2489) i nt 
2} 
~ 2 kin 3 atta [Hs Bees) 


+ ep 21(RF) Sobol vee 
Jo(R?) 


The average velocity has the following form wren 


returnedto real notation: 
_ K 2d ) 
Vitjee = 11 = == 2 t 
n | mid 5 Fer inn 


The @bove pressure history of the comoressible volume 
and the average velocity has gouping of parameters which 
appear similar to those found in the rectansular co- 


ordinate analysis. 





(1) 


(2) 


(3) 


(4) 


(5) 
(6) 


(7) 
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